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SUMMARY 


The problem of controlling large, flexible space systems has been the 
subject of considerable research. Many approaches to control system 
synthesis have been evaluated using computer simulation. In several cases, 
ground experiments have also been used to validate system performance under 
more realistic conditions. There remains a need, however, to test 
additional control laws for flexible spacecraft and to directly compare 
competing design techniques. In this paper an NASA program is discussed 
which has been initiated to make direct comparisons of control laws for, 
first, a mathematical problem, then an experimental test article is being 
assembled under the cognizance of the Spacecraft Control Branch at the NASA 
Langley Research Center with the advice and counsel of the IEEE Subcom- 
mittee on Large Space Structures. The physical apparatus will consist of a 
softly supported dynamic model of an antenna attached to the Shuttle by a 
flexible beam. The control objective will include the task of directing 
the line-of-sight of the Shuttle/antenna configuration toward a fixed 
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target, under conditions of noisy data, limited control authority and 
random disturbances. The open competition started in the early part 

of 1984. Interested researchers are provided information intended to 
facilitate the analysis and control synthesis tasks. A workshop is planned 
for early December at the NASA Langley Research Center to discuss and 
compare results. 


INTRODUCTION 

Many future spacecraft will be large and consequently quite flexible. 
As the size of antennae is increased, the frequencies of the first flex- 
ible modes will decrease and overlap the pointing system bandwidth. It 

will no longer be possible to use low gain systems with simple notch 
filters to provide the required control performance. Multiple sensors and 
actuators, and sophisticated control laws will be necessary to ensure 
stability, reliability and the pointing accuracy required for large, 
flexible spacecraft. 

Control of such spacecraft has been studied with regard given to 
modeling, order reduction, fault management, stability and dynamic system 
performance. Numerous example applications have been used to demonstrate 
specific approaches to pertinent control problems. Both computer simula- 
tions and laboratory experiment results have been offered as evidence of 
the validity of the approaches to control large, flexible spacecraft. 
Concerns remain, however, because of the chronic difficulties in control- 
ling these lightly damped large-scale systems. Because of these concerns 
and because of the desire to offer a means of comparing technical 
approaches directly, an NASA/IEEE Design Challenge is being offered. An 
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experimental teat article la being assembled under the cognisance of the 
Spacecraft Control Branch at the NASA Langley Research Center with the 
advice and counael of the IEEE ( C OLSS) Subcommittee on Urge Space 
Structures. This Spacecraft Control Laboratory Experiment (SCOLE) will 
serve as the focus of a design challenge for the purpose of comparing 
directly different approaches to control synthesis, modeling, order 
reduction, state estimation and system identification. 

The configuration of the SCOLE will represent a large antenna attached 
to the Space shuttle arbiter by a flexible beam. This configuration was 
chosen because of Its similarity to proposed space flight experiments and 
proposed space-based antenna systems. This paper will discuss the "Design 
Challenge" In terms of both a mathematical problem and a physical experi- 
mental apparatus. The SCOLE program Is not part of any flight program. 


SYMBOLS 


a 

A 

c 

d 

e 

E 

f 

F4 

g 

GI 

I 


acceleration vector ft/sec 2 
beam cross section area 
observation matrix 

noise contaminating direction cosine matrix measurements 

line-of-sight error 

modulus of elasticity 

concentrated force expressions 

force vector 


concentrated moment expressions 
torsional rigidity 


moment of inertia matrix for 


entire Shuttle/antenna configuration 


3S0 



Il 

14 
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i y 

L 

Mi 

M4 

M D 

m 

m l 

m^ 
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V 1 

v 4 

u 4> 

u 0 

U\jr 

X,Y,Z 

A 
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moment of inertia matrix, Shuttle body 

moment of inertia matrix, reflector body 

beam cross section moment of inertia, roll bending 

beam cross section moment of inertia, pitch bending 

beam polar moment of inertia, yaw torsion 

length of the reflector mast, beam 

control moment applied to the Shuttle body 

control moment applied to the reflector body 

disturbance moment applied to the Shuttle body 

mass of entire Shuttle/antenna configuration 

mass of Shuttle body 

mass of reflector body 

mass density of beam 

beam position variable 

direction cosine matrix, Shuttle body () _ * T, ()„^ , , , 

J earth 1 Shuttle body 

direction cosine matrix, reflector body () , = T, () 

J earth 4 reflector 

body 

inertial velocity, Shuttle body 
inertial velocity, reflector body 
lateral deflection of beam bending in y-z plane 
lateral deflection of beam bending in x-z plane 
angular deflection of beam twisting about z axis 
position variables 

displacement of proof-mass actuator 

line-of-sight pointing requirement 

noise contaminating angular velocity measurements 
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pitch, roll, heading 
c damping ratio 

T n ° l8e laminating acceleration measurements 

“1 angular velocity of Shuttle body 

w 4 angular velocity of reflector body 


DISCUSSION 

The objective of the NASA- IEEE Design challenge concerning the control 

flexible spacecraft la to promote direct comparison of different 

pp oaches to control, state estimation and systems Identification. The 

design challenge has principal parts, the first nslng a mathematical model 

end the second u.l„ laboratory experimental apparatus. The specific parts 

of the Spacecraft Control Laboratory Experiment (STOLE) program .111 be 
discussed in detail. 




The primary control task Is to rapidly slew or change the line-of- 
slght of an antenna attached to the space Shuttle arbiter, and to settle or 
damp the structural vibrations to the degree reared for precise pointing 
of the antenna. The objective mill be to minimise the time retired to 
slew and settle, until the antenna line-of-slght remains within the 

“ gl ' ! S ‘ * S ' COndary C<mtr ° 1 ‘o to change direction during the 

"on-target" phase to prepare for the next slew maneuver. The objective Is 

to change attitude and stabilise as quickly as possible, while keeping the 
line-of-sight error less than 6 . 
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Math Model Dynamics 


The Initial phase of the design challenge will use a mathematical 
model of the Shuttle orblter/antenna configuration. It is necessary to 
obtain a balance, of course, between complex formulations which might be 
more accurate and simplified formulations which ease the burden of 
analysis. 

The dynamics are described by a distributed parameter beam equation 
with rigid bodies, each having mass and inertia at either end. One body 
represents Space Shuttle orbiter; the other body is the antenna reflector. 
The equations for the structural dynamics and Shuttle motion are formed by 
adding to the rigid-body equations of motion, beam-bending and torsion 
equations. The boundary conditions at the ends of the beam contain the 
forces and moments of the rigid Shuttle and reflector bodies. The 
nonlinear kinetmatics couples the otherwise uncoupled beam equations. 
Additional terras represent the action of two, 2— axis proof-mass actuators 
at locations on the beam chosen by the designer. 

The rigid-body equations of motion for the Shuttle body are given by: 

"i ■ - ^vri + "i ♦ «d + v> 
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Similarly, for the reflector body. 


,-l 


"4 ‘ - : 4 <“4*4“* + M 4 * " 8 , 4 ’ 


- fil. V 

4 m. 


The direction cosine matrices defining the attitudes of the Shuttle and 
reflector bodies are given by: 


• T /"w T 

x = - CO T 

1 11 


• X X 

T 4 - - V 4 


The direction cosine matrices defining the attitudes of the Shuttle and the 
reflector bodies are related to the beam end conditions. 
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The equations of motion for the flexible beam-like truss connecting the 
reflector and Shuttle bodies consist of standard beam bending and torsion 
partial differential equations with energy dissapative terras which enable 
damped modes with constant characteristics for fixed, though dynamic, end 
conditions. The system of equations can be viewed as driven by changing 
end conditions and forces applied at the locations of the proof-mass 
actuators. 


ROLL BEAM BENDING: 
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PITCH BEAM BENDING: 
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{REFLECTOR BODY FORCE} 


{SHUTTLE BODY FORCE} 
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The angular velocity of the reflector body is related to the Shuttle body 
by: 
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The line-of-sight error described in figure 2 is affected by both the 
pointing error of the Shuttle body and the misalignment of the reflector 
due to the deflection of the beam supporting the reflector. The line-of- 
sight is defined by a ray from the feed which is reflected at the center of 
the reflector. Its direction in the Shuttle body coordinates is given by: 


-«r + r f + 2 

[ R X - v • \] 

R r _ R p ~ 2 

r a r r r f 

• r a| 


where: 


Rp is the feed location (3.75, 0, 0) 

Rr is the location of the center of the reflector (18.75, -32.5, 
-130) 

Ra is a unit vector in the direction of the reflector axis in 
Shuttle body coordinates 


The vector R A can be related to the direction cosine attitude matrices 
for the Shuttle body, Tj , and the reflector body, T^ , by 



The relative alignment of the reflector to the Shuttle body is given by 
T 

T^T^ which is a function of the structural deformations of the beam. 
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The line-of-sight error, e, is the angular difference between the 
target direction, given by the unit vector, Df» and the line - of - sight 
direction in Earth axes, TiR L os* 

e - ARCSIN |D T X or ARCSIN (d^R^ 

Computer programs are available which generate time histories of the 
rigid body and the mode shapes and frequencies for the body-beam-body 
configuration for "pitch” bending, "roll" bending and "yaw" twisting. 
Since the modes are based on solving explicitly the distributed parameter 
equations (without damping and without kinematic coupling) there is no 
limit to the number of modal characteristic sets that can be generated by 
the program. It will be the analyst's decision as to how many modes need 
to be considered. 


Laboratory Experiment Description 

The second part of the design challenge Is to validate in the 

laboratory, the system performance of the more promising control system 
designs of the first part. The experimental apparatus will consist of a 
dynamic model of the Space Shuttle orbiter with a large antenna reflector 
attached by means of a flexible beam. The dynamic model will be exten- 

sively instrumented and will have attached force and moment generating 
devices for control and for disturbance generation. A single, flexible 
tether will be used to suspend the dynamic model, allowing complete angular 
freedom in yaw, and limited freedom in pitch and roll. An inverted 

position will be used to let the reflector mast to hang so that gravity 

effects on mast bending will be minimized. The dynamics of the laboratory 
model will of necessity be different from the mathematical model discussed 
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earlier. 


Design Challenge, Part One 


For part one of the design challenge, the following mathematical 
problem is addressed. Given the dynamic equations of the Shuttle/antenna 
configuration, what control policy minimizes the time to slew to a target 
and to stabilize so that the line-of-sight (LOS) error is held, for a time, 
within a specified amount, 6. During the time that the LOS error is 
within 6, the attitude must change 90° to prepare for the next slew 
maneuver. This was previously referred to as the sescondary control task. 
The maximum moment and force generating capability will be limited. Advan- 
tage may be taken of selecting the most suitable initial alignment of the 
Shuttle/antenna about its assigned initial RF axis, line-of-sight. 

Random, broad band-pass disturbances will be applied to the configuration. 
Two proof-mass, force actuators may be positioned anywhere along the beam. 
The design guidelines are summarized below: 

1. The initial line-of-sight error is 20 degrees. 

e(o) =* 20 degrees 

2. The initial target direction is straight down. 



3. The initial alignment about the line-of-sight is free to be chosen 
by the designer. Advantage may be taken of the low value of 

moment of inertia in roll. The Shuttle/antenna is at rest 
initially, 

4. The objective is to point the line-of-sight of the antenna and 
stabilize to within 0.02 degree of the target as quickly as 
possible. 

5 = 0.02 degree 

3G0 
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5. Control moments can be applied at 100 Hz sampling rate to both the 
Shuttle and reflector bodies of 10,000 ft-lb for each axis. The 

commanded moment for each axis Is limited to 10,000 ft-lb. The 

actual control moment's response to the commanded value is 
first-order with a time constant of 0.1 second. 

For the rolling moment applied to the Shuttle body: 


- 10 < 


- Yi, 


. < 10 
command — 


- 0.1 


- 0 . 1 , 


M x,l <n + 0 * e Yl 00 + U ' e ' > Yl,comm«>d (n) 


Equations for other axes and for the reflector body are similar. 

6. Control forces can be applied at the center of the reflector in 
the X and Y directions only. The commanded force in a 

particular direction is limited to 800 lbs. The actual control 
force's response to the commanded value is first-order with a 
response time of 0.1 second. 

For the side for applied to the reflector body: 


-800 < F v . < 800 

— Y , command — 


F (n + 1) = e' 0 * 1 F (n) + (1 - e °* 1 ) F .(n) 

Y Y Y , command 


Equations for X-axis are similar. 

7. Control forces using two proof-mass actuators (each having both 
X and Y axes) can be applied at two points on the beam. The 

strokes are limited to ± 1 ft, and the masses weight 10 lbs each. 
The actual stroke follows a first-order response to limited 
commanded values. 


jet 



For the X axis of the proof- mass actuator at 


s 2 : 


-1 < A <1 

— X , 2 , command — 


- 0.1 


V ( " + *> ' 4 x,2<"> * U - e-°-‘) \ i2 ' Comland M 


Equations for other axes and locations are similar. 

8. The inertial attitude direct ton cosine matrix for the Shuttle body 
lags in time the actual values by 0.01 second and are made at a 
rate of 100 samples per second. Each element of the direction 
cosine measurement raatarix is contaminated by additive, 
uncorrelated Gaussian noise having an rras value of 0.001. The 
noise has zero mean. 



d u<"> 

d 12 (n) 

d l 3 (n) 

T ,(n + 1) = T (n) + 

s, measured s,true v ' 

d 21 (n) 

d 22 ( "> 

d 23^ n > 

where: 

d 31<"> 

d 32 (n) 

d 33 (n) 

E{d i:j (n)} = 0 

E{d (n)d kL (n)} = 0 

for i 

* k or 

j * L 

E{d (n)d (n + k)} = 0 

for k 

* 0 


= 1 .001 j 2 

for k 

= 0 
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9. The angular velocity measurements for both the Shuttle and 
reflector bodies pass through a first-order filter with 0.05 sec 
time constant and lag in time the actual values by 0.01 second and 
are made at a rate of 100 samples per second. Each rate 
measurement is contaminated by additive, Gaussian, uncorrelated 
noise having an rms value of 0.02 degree per second. The noise 
has zero mean. 

For example: 

W 1 ,X, measured^" + ^ “l ,X, f iltered* n) + E l,X (n) 

E{ e i ^(n) Gj ^(n + k)} = 0 for k * 0 

«* ( .02) 2 for k = 0 

where 

GU = - 20 Cl) + 20 cu 

1,X, filtered 1,X, filtered l,X,true 

10. Three-axis accelerometers are located on the Shuttle body at the 
base of the mast and on the reflector body at its center. Two- 

axes (X and Y) accelerometers are located at intervals of 
10 feet along the mast. The acceleration measurements pass 

through a first-order filter with a 0.05 second time constant and 
lag in time the actual values by 0.01 second, and are made at a 
rate of 100 samples per second. Each measurement is contaminated 
by Gaussian additive, uncorrelated noise having an rms value of 
0.05 ft/sec^. 
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For example : 


a l ,x, measured^" + ^ *1 ,X, filtered^ + T 1,X^ 

e{t (n) T (n + k)} = 0 for k * 0 

i )X i 

= ( ,05) 2 for k = 0 

where: 

3# 1 ,X, filtered = " 20 3 1 ,x , filtered + 20 ^l.X.true 

11. Gaussian, uncorrelated step-like disturbances are applied 
100 times per second to the Shuttle body in the form of 3-axes 
moments, having rms values of 100 ft-lbs. These disturbances 

have zero mean. 

For example: 

E { M D x (n) M d x (n + k)} = 0 for k * 0 

= (100) 2 for k = 0 

In summary, the designer's task for part one is to: (1) derive a 

control law for slewing and stabilization, coded in FORTRAN; (2) select an 
initial attitude in preparation for slewing 20 degrees; and (3) select two 
positions for the 2-axes proof-mass actuators. An official system 
performance assessment computer program will be used to establish the time 
required to slew and stabilize the Shuttle/antenna configuration. 
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Design Challenge, Part Ttoo 


As in part one, the task is to minimize the time to slew and stabilize 
a Shuttle/antenna configuration. The difference is that in part two of the 
design challenge, a physical laboratory model will be used instead of the 
dynamic equations of part one. The constraints on total moment and force 

generation capability will apply to part two, as for part one. Again, the 

analyst may select the initial alignment about the assigned initial RF 

line-of-sight. Disturbances will be injected into the Shuttle/antenna 
model. The designer’s task will be similar to that for part one. 

CONCLUDING REMARKS 

A Design Challenge, in two parts, has been offered for the purpose of 

comparing directly different approach to controlling a flexible 

Shuttle/antenna configuration. The first part of the design challenge uses 
only mathematical equations of the vehicle dynamics; the second part uses a 
physical laboratory model of the same configuration. The Spacecraft 
Control Laboratory Experiment (SCOLE) program is being conducted under the 
cognizance of the Spacecraft Control Branch at the NASA Langley Research 
Center. The NASA/ IEEE Design Challenge has the advice and counsel of the 
IEEE-COLSS Subcommittee on Large Space Structures. Workshops will be held 
enable investigators to compare results of their research. 
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The moment of inertia becomes: 
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CONTROL MOMENTS, DISTURBANCES 





AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE 



Roll bending mode $ 1. Frequency^ .32 Hz 



Roll bending mode # 2. Frequency= 1.29 Hz 



Roll bending mode § 3. Frequency= 480 Hz 



Roll bending mode # 4 Frequency= 12.29 Hz 



Roll bending mode § 5. Frequency= 23.68 Hz 



Roll bending mode § 6. Frequency 38.89 Hz 



Roll bending mode # 7. Frequency= 57.90 Hz 



Roll bending mode # 8. Frequency^ 80.72 Hz 


Figure 4.i.- plots of normalized roll bonding mode shapes 
for SCOLE configuration. 
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wiPLiTuDE anPLiTuoe amplitude amplitude amplitude 



Pitch bending mode # 1. Frequency= .29 Hz 



Pitch bending mode # 2. Frequency= 1.65 Hz 



Pitch bending mode # 3. Frequency= 4.97 Hz 



Pitch bending mode # 4 Frequency^ 12.36 Hz 



Pitch bending mode # 5. Frequency= 23.72 Hz 



Pitch bending mode # 6. Frequency= 38.91 Hz 



Pitch bending mode # 7. Frequency= 57.92 Hz 



Pitch bending mode # 8. Frequency= 80.73 Hz 


Figure 4b,- Plots of normalized pitch bending mode shapes 
for SCOLE configuration. 
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Torsional mode # 1. Frequency- .53 Hi 
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Torsional mode § 2. Frequency= 45.12 Hz 



Torsional mode # 4. Frequency= 135.35 Hz 



Torsional mode # 6. Frequency = 225.57 Hz 
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Torsional mode # 7, 
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Frequency = 270.69 Hz 



_ \y 

Torsional mode # 8, 
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Frequency= 315.80 Hz 


Figure 4c,- Plots of normalized torsional mode shapes for 
SCOLE configuration . 
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